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Abstract. It is generally accepted that scale-free networks is prone to epidemic spreading allowing the 
onset of large epidemics whatever the spreading rate of the infection. In the paper, we show that disease 
propagation may be suppressed in particular fractal scale-free networks. We first study analytically the 
topological characteristics of a network model and show that it is simultaneously scale-free, highly clus- 
tered, "large- world" , fractal and disassortative. Any previous model does not have all the properties as 
the one under consideration. Then, by using the renormalization group technique we analyze the dynamic 
susceptible-infected-removed (SIR) model for spreading of infections. Interestingly, we find the existence of 
an epidemic threshold, as compared to the usual epidemic behavior without a finite threshold in uncorre- 
cted scale- free networks. This phenomenon indicates that degree distribution of scale- free networks does 
not suffice to characterize the epidemic dynamics on top of them. Our results may shed light in the un- 
derstanding of the epidemics and other spreading phenomena on real-life networks with similar structural 
features as the considered model. 



PACS. 89. 75. He Networks and genealogical trees 
Stability and fragmentation of clusters 



1 Introduction 

' In the past ten years, there has been a considerable in- 
terest in characterizing and understanding the topological 
properties of networked systems 1,2,3,4,5 . It has been es- 
tablished that scale- free behavior [Hj is one of the most fun- 
damental concepts constituting our basic understanding of 
the organization of many real- world systems in nature and 
society. This scale-free property has a profound effect on 
almost every aspect on dynamic processes taking place on 
networks, including robustness [7], percolation [ElE], syn- 
chronization [TT], games [T2], epidemic spreading [13 ,14, 
ITS] , and so on. For instance, for a wide range of scale- free 
networks, there is no existence of an epidemic threshold, 
even infections with low spreading rate will prevail over 
the entire population in these networks p31[l4j[l5] • This 
radically changes the conclusion drawn from classic dis- 
ease modeling [16] , 

Recently, it has been discovered that many real-life 
networks, such as the WWW, metabolic networks, and 
yeast protein interaction networks have self-similar prop- 
erties and exhibit fractal scaling [P? | I18 | I19 | I20|. The frac- 
tal topology can be characterized through two exponents: 
fractal dimension ds and degree exponent of the boxes 
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dk , which can be obtained by box-counting algorithm |17[ 
18,21,22 . The scaling of the minimum number of boxes 
Nb of linear size £b needed to cover the network with 
node number N defines the fractal dimension ds, namely 
Nb/N ~ . Analogously, the degree exponent of the 

£g dk , where 
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boxes dk is identified through fes(^s)/fc/iu6 
fc.B(^.B) is the number of outgoing links from the box as a 
whole, and khub t ne largest node degree inside the box. 
Fractal networks are all self-similar, which means that 
fractal scale-free networks present the property of scale- 
invariance of degree distribution, P(k) ~ k i.e., the 
exponent 7 remains the same for different box sizes |17) . 
In self-similar scale-free networks, the three indexes 7, ds 
and dk satisfy the following relation: 7 = 1 + ds/dk [17] . 

As a fundamental property, topological fractality re- 
lates to many respects of network structure and function. 
Recent authors have shown that the correlation between 
degree and betweenness centrality of nodes is much weaker 
in fractal network models in comparison with non-fractal 
models [23]- It has been also shown [18lfT9ir24] that frac- 
tal scale-free networks are not assortative, this disassor- 
tativity feature together with fractality makes such scale- 
free networks more robust against intentional attacks on 
hub nodes, as compared to the very vulnerable non-fractal 
scale- free networks [T5]. In addition to the distinction in 
the robustness, fractal networks exhibit poorer synchro- 
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nizability than non-fractal counterparts [24] • Although a 
lot of jobs have been devoted to characterizing fractal net- 
works 25,26,27,28,29,30], it is still of current interest to 
model fractal topology and seek a better understanding of 
its consequences on dynamic processes. 

In this paper, we relate fractality to dynamics of dis- 
ease spread in deterministic networks. Deterministic graphs 
have strong advantages. For example, they allow to com- 
pute analytically their properties, which have played a sig- 
nificant role, both in terms of explicit results and a guide 
to and a test of simulated and approximate methods [314 
13211331 13 ^ 131 ^ ^ ^ 

We first introduce a deterministic family of fractal graphs. 
From the viewpoint of complex networks, we determine 
accurately the topological characteristics of a particular 
graph and show that it is simultaneously scale- free, highly 
clustered, fractal and disassortative, but lacks the small- 
world property. Then we define a dynamic susceptible- 
infected-removed (SIR) model j!6j on the two dimensional 
fractal graph to study the effect of fractality on disease 
spreading. By mapping the SIR model to a bond percola- 
tion problem, we found that there is an existence of finite 
epidemic threshold. Thus, the transmission rate needs to 
exceed a critical value for the disease to spread and pre- 
vail, which shows that the fractal networks are robust to 
infection. 



2 Network construction and topologies 

This section is devoted to the construction and the rel- 
evant structural properties of the studied network, such 
as degree distribution, clustering coefficient, average path 
length (APL), fractality, and correlations. 




Y Y Y 



q=l q=2 

Fig. 1. (Color online) Iterative construction method of the 
fractal networks. Each iterative link is replaced by a connected 
cluster on the right-hand side of the arrow. The blue links are 
noniterated ones. 




t=2 t=3 



Fig. 2. (Color online) Scheme of the growth of the network 
for the particular case of q = 1. Only the first three iterative 
processes are shown. 



are generated as follows: For t = 0, Fo. q is an iterative 
edge connecting two nodes. For t > 1, F t . q is obtained from 
Ft~i,q- We replace each existing iterative edge in i*t-i,g by 
a connected cluster of edges on the right of Fig. [TJ where q 
denotes the number of the newly-created noniterated links 
induced by the iterative edge. The growing process is re- 
peated t times, with the family of fractal graphs obtained 
in the limit t — > oo. Figure [2] shows the growing process of 
the network for the particular case of q = 1. Note that all 
generated networks have qualitatively similar properties. 
In what follows we focus on the special case of q — 1 and 
denote it as F t after t iterations. 

Next we compute the numbers of total nodes (vertices) 
and links (edges) in Ft- Notice that there are two types 
of links (i.e., iterative links and noniterated links) in the 
network. Let L v (t), Li(t) and L n (t) be the number of new 
vertices, iterative links, and noniterated links created at 
step t, respectively. Since all old iterative links are not 
preserved in the growing process, thus Li(t) is in fact the 
total number of iterative links at time t. Note that each of 
the existing iterative links yields two nodes connected by 
one noniterated link, and the addition of each new node 
leads to two iterative links. By construction, for t > 1, we 
have 

Li{t)=ALi{t-l), 
L v (t) = 2Li(t-l), (1) 
L n (t)=Li(t-l). 

Considering the initial condition L V (Q) — 2, Li(0) = 1, 
and L n (Q) = 0, it follows that 

L v (t) = 2-4*-\ 

Li(t) = 4*, (2) 



Thus the number of total nodes N t and edges E t present 
2.1 Construction algorithm a t s t e p i } s 

The proposed class of fractal networks is constructed in t , 

an iterative way as shown in Fig.[TJ Let F t ^ (t > 0, q > 1) Nt= L v (ti) = - ' - — - (3) 

denote the networks after t iterations. Then the networks t ._ 3 
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and 



2.3 Clustering Coefficient 



4 t+1 - 1 

E t =Li{t) + J^L n (ti) = , (4) 

respectively. The average degree after t iterations is 

_2_E 1 _ 2(4'+* - 1) 
Wi " N t ~ 2-4' + 4 ' (5) 

which approaches 4 in the infinite t limit. 
2.2 Degree distribution 

When a new node u is added to the network at step t u 
(t u > 1 ) , it has three links, among which two are iterative 
links and one is noniterated link. Let Li(u, t) be the num- 
ber of iterative links at step t that will create new nodes 
connected to the node u at step t + 1. Then at step t u , 
Li(u, t u ) = 2. From the iterative generation process of the 
network, one can see that at any subsequent step each it- 
erative link of u is broken and generates two new iterative 
links connected to u. We define k u (t) as the degree of node 
u at time t, then the relation between k u (t) and Li(u,t) 
satisfies: 

k u (t)=Li(u,t) + l, (6) 

where the last term 1 represents the only noniterated link 
of node u. Now we compute Li(u,t). By construction, 
Li(u,t) = 2Li(u,t — 1). Considering the initial condition 
Li(u,t u ) — 2, we can derive Li(u,t) — 2 t ~ tu+1 . Then at 
time t, the degree of vertex u becomes 



k u (t) = 2 



t-t u +i 



1. 



(7) 



It should be mentioned that the initial two nodes cre- 
ated at step have a little different evolution process from 
other ones. Since the initial two nodes have no noniterated 
link, we can easily obtain that at step t, for either of the 
initial two nodes, its degree just equals the number of it- 
erative links connecting it, both of which are 2*. 

Equation ([7]) shows that the degree spectrum of the 
network is discrete. It follows that the cumulative degree 
distribution [3] is given by 



^cum(^) — 



E 

T<t u 



L v (t) _ 2-4*" +4 



Nt 



2-4* 



(8) 



Substituting for t u in this expression using t,, 

ln(fc-l) • 

4^ S lves 



t + 1 



-fcum(^) 



2-4 4 -4(fc-l)-( ln4 / ln2 )+4 



2 • 4* + 4 

When t is large enough, one can obtain 

Pcum(£0=4(fc-1)- 2 . 



(9) 



(10) 



So the degree distribution follows a power law form with 
the exponent 7 = 3. The same degree exponent has been 
obtained in the famous Barabasi- Albert (BA) model [5] 
as well as some other deterministic models [28,24,29,30, 
I5ni52ll53l54l. 



The clustering coefficient [55] of a node u with degree k u 
is given by C u = 2e u /[k u (k u — 1)], where e u is the num- 
ber of existing links among the k u neighbors. Using the 
construction rules, it is straightforward to calculate ana- 
lytically the clustering coefficient C(k) for a single node 
with degree k. For the initial two nodes born at step 0, 
their degree is k — 2*, and the existing links among these 
neighbors is 4, all of which are noniterated links. For those 
nodes created at step cj) (0 < < i), there are only 
links that actually exist among the neighbor nodes. Fi- 
nally, for the smallest nodes created at step t, each has a 
degree of k = 3, the existing number of links between the 
neighbors of each is 2. Thus, there is a one-to-one corre- 
spondence between the clustering coefficient C(k) of the 
node and its degree k: 



C(k) = 




1) for k = 2* 

for k = 2 m + 1(2 < m < t) 
for k = 2 1 + 1 



(11) 



which is inversely proportional to k in the limit of large 
k. The scaling of C(k) ~ fc" 1 has been observed in many 
real- world scale- free networks |39j . 

Using Eq. (fTTj) . we can obtain the clustering C t of 
whole the network at step t, which is defined as the aver- 
age clustering coefficient of all individual nodes. Then we 
have 



MQ) 

D -l 



t-i 

E 

r=l 



(r) 2L v (t) 



Dt 



(12) 



where D r is the degree of a node at time t, which was 
created at step r, see Eq. ([7]). In the infinite network order 
limit (Nt — > oo), Eq. (| 1 2[) converges to a nonzero value 
C = 0.5435. Therefore, the average clustering coefficient 
of the network is very high. 



2.4 Fractal dimension 

As a matter of fact, the fractal graph grows as a inverse 
renormalization procedure, see Fig. [3 in reverse order. To 
find the fractal dimension, we follow the mathematical 
framework presented in Ref. [18] . By construction, in the 
infinite t limit, the different quantities grow as: 



JV t =i4JV t _i, 
fe u (t)~2fe u (t-l), 
D i = 2%_ 1)) 



(13) 



where the third equation describes the change of the diam- 
eter D t of the graph F t , where Bt is defined as the longest 
shortest path between all pairs of nodes in F t . 

From the relations provided by Eq. (fl3|) . it is clear 
that the quantities Nt, k u (t) and D t increase by a factor 
of Jn = 4, fk = 2 and fo — 2, respectively. Then between 
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any two times t\, t<i (ti < t-i), we can easily obtain the 
following relation: 



D* 2 = 2 t2 - tl ©t 1 , 
N t2 =^-^N tl , 
k u {t2) = K(t x ). 



(14) 



From Eq. (|14| . we can derive the scaling exponents in 
terms of the microscopic parameters: the fractal dimension 
is ds — l \l/fl — 2, and the degree exponent of boxes is 
dk = = 1. The exponent of the degree distribution 

satisfies 7 = 1+^ = 3, giving the same 7 as that obtained 
in the direct calculation of the degree distribution, see 
Eq. (HQ}. 

Note that in a class of deterministic models called 
pseudo- fractals, although the number of their nodes in- 
creases exponentially, the additive growth of the diame- 
ter with time implies that the networks are small world. 
These models do not capture the fractal topology found 
in diverse complex networks pllMllMll35 l l3^ 1 l ^ l4 ^ l4T | l42 l 



2.5 Degree correlation 

Degree correlation [5^[S71[5^[5^1ISnilr5llS2l lB3] is a particu- 
larly interesting subject in the field of complex networks, 
because it can give rise to some interesting network struc- 
ture effects. Degree correlation in a network can be mea- 
sured by means of the quantity, called average nearest- 
neighbor degree (ANND) and denoted as k nn (k), which 
is a function of node degree, and is more convenient and 
practical in characterizing degree correlation. ANND is 
defined by [S7] 



k nn (k)=Y J k'P{k'\k), 



(15) 



where P(k'\k) is the probability that a link from a node 
of degree k points to a node of degree k' . 

For the fractal graph considered here, one can exactly 
calculate k nn (k) . By construction, all neighbors of the ini- 
tial two nodes have the same degree 3, while for each other 
nodes with degree greater than 3, only one of its neighbor 
has the same degree as itself, all the rest neighbors have 
a degree 3. Then we have 

f k nn (k) =3 for k = 2* 

I k nn (k) = 4 - | for fc = 2 m + l(m = 2,3 •••£). ^ 

For those nodes with degree 3, it is easily to obtain 

fcnn(3) = 



Degree correlation can also be described by a Pearson 
correlation coefficient r of degrees at either end of a link. 
It is defined as [59l[60ll42ll64l 



2-V t ) 2 +T, T r=l[Lv(T)k(T,t)[k(T,t)-l)] 

3L„(t) 



4 5 4 1 
3* + 3 ~ 3 ' 2*' 



(17) 



where k(r, t) is the degree of a node at time t that was born 
at step t. Thus /c nn (3) grows linearly with time for large 
t. Eqs. (|16| and (fl7| show the network is disassortative. 



(fc)(fc 2 fc nn (fc))-(fc 2 ) ; 
(fc)(fc 3 ) - (fc 2 ) 2 



(18) 



We can easily see that for t > 1, r of Ft is always negative, 
indicating disassortativity. 



2.6 Average path length 

We represent all the shortest path lengths of F t as a matrix 
in which the entry dij is the shortest path from node i 
to j. A measure of the typical separation between two 
nodes in F t is given by the average path length ( APL) dt , 
also known as characteristic path length [5] , defined as the 
mean of geodesic lengths over all couples of nodes. APL is 
relevant in many fields regarding real- life networks and has 
received much attention [65.66,67,168). In the Appendix, 
we have obtained exact analytic expression for dt, which 
reads 



d t = 



(16 • 2* + 21)16* + (21* - 27)8* + 75 • 4* + 119 • 2* - 15 



21(2 + 5 -4* + 2 -16') 



(19) 



For large t, d t — ► Jj- • 2*. Note that in the infinite t limit, 

_ 1/2 

N t ~ 4*, so the APL scales as dt ~ N t , which indicates 
that the network is not a small world. 

Thus, we have shown that d t has the power-law scal- 
ing behavior of the number of nodes Nt, which is simi- 
lar to that of two-dimensional regular lattice [69]. This 
phenomenon is not hard to understand. Let us look at 
the scheme of the network growth. Each next step in the 
growth of F t doubles the APL between a fixed pair of 
nodes (except those small number of pairs directly con- 
nected by a noniterated link), while the total number of 
nodes increases four-fold (asymptotically, in the infinite 
limit of t), see Eq. ([3]). Thus the APL d t of F t grows as a 
square power of the node number in the network. 



3 SIR model on the network 

As discussed in previous section, the network exhibits many 
interesting properties, i.e., it is simultaneously scale-free, 
highly clustered, "large-world" , fractal and disassortative, 
which is not observed in uncorrelated networks with purely 
random wiring. Therefore, it is worthwhile to investigate 
the processes taking place upon the model to find the dif- 
ferent impact on dynamic precesses compared with uncor- 
related networks. In what follows we will study the SIR 
model of epidemics, which is one of the first issues to be 
explored in the literature on complex networks [13, 141ll5j. 

In the standard SIR model [16], each node of the net- 
work represents an individual and each link is the con- 
nection along which the individuals interact and the epi- 
demic can be transmitted. This model describes diseases 
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resulting in the immunization or death of infected indi- 
viduals, and assumes that each individual can be in one 
of three possible states, namely, susceptible, infected, and 
removed. The disease transmission on the network is de- 
scribed in an effective way: At each time step, each suscep- 
tible node is infected with probability A, if it is connected 
to one or more infected nodes; at the same time, each 
infected individual becomes removed with probability 1, 
therefore it can not catch the infection again. 

The SIR model is equivalent to a bond percolation 
problem with bond occupation probability A 70, 71\ . More- 
over, the size of the outbreak is just the size of the giant 
component. In our case the percolation problem can be 
solved using the real-space renormalization group tech- 
nique 72,73,74,30 , giving exact solution for the interest- 
ing quantity of epidemic threshold. Let us describe the 
procedure in application to the considered network. As- 
suming that the network growth stops at a time step 
t — > oo, when the network is spoiled in the following way: 
for a link present in the undamaged network, with the 
probability A we retain it in the damaged network. Then 
we invert the transformation in Fig. 2 and define n = t — r 
for this inverted transformation, which is actually a dec- 
imation procedure |74j . Further, we introduce the proba- 
bility A„ that if two nodes are connected in the undam- 
aged network at r = t — n, then at the nth step of the 
decimation for the damaged network, there exists a path 
between these vertices. Here, Ao = A. We can easily obtain 
the following recursion relation for X n 

A n+1 = X 5 n + 5A*(1 - A„) + 8X 3 n (l - A„) 2 + 2X 2 n (1 - A„) 3 
= 2A 5 „ - 5\ 4 n + 2\l + 2A 2 . (20) 

Equation (|20[) has an unstable fixed point at A c = h, and 
two stable fixed points at A = and A = 1. 

Thus, for the SIR model the epidemic prevalence un- 
dergoes a phase transition at a finite threshold A c of the 
transmission probability. If infection rate A > A c , the dis- 
ease spreads and infects a finite fraction of the popula- 
tion. On the other hand, when A < A c , the total number 
of infected individuals is infinitesimally small in the limit 
of very large populations. The existence of an epidemic 
threshold in the present network is compared to the re- 
sult for uncorrelated scale-free networks, where arbitrarily 
small infection rate A shows finite prevalence [T3Ill"4"Ill"5] . 

Why uncorrelated scale- free networks are prone to epi- 
demics spreading, while the present fractal disassortative 
scale-free network can suppress disease propagation? This 
may be explained as follows. In non-fractal uncorrelated 
scale-free topologies, the hubs are connected and form a 
central compact core, such that the infection of a few of 
the largest hubs has catastrophic consequences for the net- 
work. For the fractal network, it is disassortative and self- 
similar, which do not allow the presence of direct con- 
nections between hubs, i.e., hubs are more dispersed (see 
Fig. Thus, the mixture of fractal and disassortative 
properties significantly provide protection against disease 
spreading. This could provide insight into explaining why 
some real-life networks have evolved into a fractal and dis- 
assortative architecture [T71IT8]. 



4 Conclusions and discussion 



To conclude, we have investigated a class of deterministic 
graphs from the viewpoint of complex networks. The de- 
terministic self-similar construction allow us to derive an- 
alytic exact expressions for the relevant features. We have 
shown that the graphs simultaneously exhibit many in- 
teresting structural characteristics: power-law degree dis- 
tribution, large clustering coefficient, 'large-world' phe- 
nomenon, fractal similar structure, negative degree cor- 
relations. The simultaneous existence of scale-free, high 
clustering, and 'large-world' behaviors is compared with 
previous network models. 

Moreover, we have studied the SIR model in the graph 
under consideration. We have presented the presence of a 
finite epidemic threshold in the finite network size limit, 
showing that being prone to disease spreading is not an 
intrinsic property of scale-free networks. The ability of 
suppressing epidemic spreading may be owing to its in- 
herent topologies. Thus, our research may be helpful for 
designing real networks resistent to epidemic outbreaks, 
and for the better understanding of the role that network 
structure plays in the spread of disease. 

Although we have studied only a particular network 
corresponding to q = 1, in a similar way, one can easily 
investigate other networks (i.e., q > 2 cases) with vari- 
ous values of 7 and ds , and their general properties such 
as 'large-world' behavior, high clustering coefficient, and 
disassortative phenomenon are similar. Analogously, one 
can also analyze the SIR model on top of these networks. 
There is an existence of a different finite epidemic thresh- 
olds for all the cases, which depend on the parameter q. 
We speculate that the fractal property of the graphs deter- 
mines the presence of the threshold of disease transmission 
for SIR model that can be mapped to a bond percolation 
as in the <1b = '"[^ dimensional regular lattices In 
the end, we should mention that since most real-world net- 
works are stochastic, it would be interesting to construct 
random network models displaying similar structural fea- 
tures as the deterministic model studied in the present 
work. 
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The last two terms at the end compensate for the over- 
counting of certain paths: the shortest path between B 
and D, with length 2 t+1 , is included in A]' 2 and A 3A ; 
the shortest path between A and C, with unit length 1, 
is included in both A\ ,A and A 2 ' 3 . 

By symmetry, A}; 2 = A^ 4 - 
A 2,3 , so that 



A 1,3 



A 



and A\ ,A 



A t = 2A^ + 2Z\ t 1 ' 3 + 2A\ 



->t+i 



(26) 



Fig. 3. (Color online) Second construction method of the net- 
work that highlights self-similarity: The graph after t + 1 con- 
struction steps, Ft+i, is composed of four copies of Ft denoted 
as F" (a = 1,2,3,4), which are connected to one another as 
above. The blue dashed link is noniterated edge. 



Appendix A: 
length 



Derivation of the average path 



By definition, the APL for F t is defined as follows 

A 



N t (N t - l)/2 ' 



(21) 



(22) 



where 

D t = di 3 
ieF t ,jeF t ,i^j 

denotes the sum of the chemical distances between two 
nodes over all pairs, and dij is the chemical distance be- 
tween nodes i and j. The network has a self-similar struc- 
ture allowing one to calculate d t analytically. As shown 
in Fig. [3j the lattice F t +i may be obtained by the jux- 
taposition of four copies of Ft, which are labeled as Lf, 
a = 1,2,3,4. Then we can write the sum A+i as 



D 



t+i 



4 A 



(23) 



where A t is the sum over all shortest paths whose end- 
points are not in the same F t branch. The solution of 
Eq. (21 is 



t-i 



A = 4'- 1 A 



(24) 



The paths that contribute to A t must all go through at 
least one of the four edge nodes (A, B, C, D, see Fig. [3]) 
at which the different F t branches are connected. The an- 
alytical expression for A tl called the crossing paths, is 
found below. 

Denote A^ as the sum of all shortest paths with end- 
points in F" and Ff. If F" and Ff meet at an edge node, 
A"'^ 3 rules out the paths where either endpoint is that 
shared edge node. If F" and Ff do not meet, Af' 13 ex- 
cludes the paths where either endpoint is any edge node. 
Then the total sum A t is 



A t =A 



1.2 



,1,3 



\1,4 



A 



2.3 



\2A 



A 3 ' 4 - 2 



t+i 



1, 

(25) 



where A]' 2 is given by the sum 
Al' 2 = d v 

hj^A 

= ^ ( d iA + d A j) 

ieFl, jeF? 
hj^A 

= (N t - 1) ^ d 'A + {Nt - 1) dA 3 

= 2(N t d i A > ( 2? ) 

where J2ieF 1 d iA = J2j£F 2 d Aj has been used. To find 
^2 i£F i diA, we examine the structure of the graph at the 
ith level. In F^, there are v t (m) points with diA = 
where 1 < m < 2 , and vt(ra) can be written recursively 
as 



u t (m) 



2' if m is odd , 

^t-i(ir) ifm is even. 



(28) 



We can write ^2 i€F i diA in terms of ^t(m) as 



2' 



iEF} 



(29) 



Eqs. (f28|) and (|29[) relate bt and bt-i, which allow one to 
resolve bt by induction as follow: 



b t = ^(2fc-l)-2 t + ^2fc-^ t _i(fc) 



k=l 



k=l 



2 3t ~ 2 + 2b t -i = -2*(2 + 4*), 
3 



(30) 



where b 1 = t^i(l) + 2^(2) = 4 has been used. Substituting 
§( 



Eq. (J30J) and N t = |(2 + 4*) into Eq. (27]), we obtain 



A]' 2 = i 2 t+1 (2 1+2t + 1) (4* + 2) 



(31) 
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Continue analogously, 
A 1,4 = E dt 3 

ieFt,i^A,D 
]£F?,j^C,D 



relation 



E (* 

d zA +d jC <2 t 



iA 



d Jc) + E dAC 
ieF^i^A 

d iA +d jC <2 t 

2J (d iD +d jD ) 

i£F£,i^A 
d iC +d jC >2 t 

E + E d * A - 2 • 2 * - 1 ' ( 32 ) 



E 



d 4 c+<i/c=2 t 



The first term equals the fourth one and is denoted by 
gt- The second and third terms are denoted by yt and 
ht, respectively. The fifth and sixth terms are equal to 
each other, both of which equal to bt + N t — 1 . The last 
two terms at the end compensate for the redundancy or 
overcounting of certain paths: the shortest paths of dco 
in X^ieF 1 di c an d d,AD in X^eF 4 d i A are redundant, both 
of which are 2*; the shortest path between dAC with unit 
length 1 is included in both XieF 1 d iC an d XjeF 4 djA- 
Therefore A\ A = 2g t + h t + y t + 2b t - 2 t+1 + 2N t - 3. One 
can compute the quantity gt as 

2*-2 2*-l-m 

9t = E E vt(m)vt(m')(m + m') 

m—l m' — l 

2 t ~ 1 -2 2 t ~ 1 -l-k 

= E E ^-i(fcK-i(fc')(2fc + 2fc') 
fc=i fc'=i 

2 t-l_ 1 2 *- 1 -fc 

+ 53 51 v t _ 1 (k)2 t {2k + 2k' - 1) 
fe=i fc'=i 

2 *-i_i 2*~ 1 -/c 

+ 53 53 2*j/ t _ 1 (fc')(2fc-l + 2fc / ) 



(33) 



fe=i fe'=i 

2 t-l_ 1 2 *-l_ fc 

+ 53 E 2 2t (2fc-l + 2fc' -1), 
fc=i fc'=i 

where the fourth term can be summed directly, yielding 
1 



-2" 2+3t (-2 + 2*)(-l + 2 t ). 



(34) 



In Eq. (|3"3"|) , the second and third terms are equal to each 
other and can be simplified by first summing over k', yield- 
ing 



2 l ^ (4* -fc 2 



(35) 



fe=r 



For use in Eq. (35), ^Li 1 ^t-l( fe ) = ^t-l ~ 2 - 0n thc 
other hand, it is easy to derive the following recursive 



53 #Vi(*) = E( 2z - i) 22 ^ 1 + E( 2z ) 2 ^*-2( fc )' 

2 = 1 

(36) 



fe=i 



2 = 1 



using which we have 



2*- L -i 



fc=i 



53 fc 2 ^ 1 (A : ) = 53^_ 1 (fc)-2 2t - 2 

.2 2t - 3 (4*-3i+17)-2 2t - 2 . (37) 
With these results, Eq. lj3"5") becomes 



fc=i 
1 

9 1 



-i+t 



-11 + 2 



1+2* 



3i) 



With Eqs. ([31) and f38), Eq. (03) becomes 



(38) 



gt =2g t - 1 + — (5 • 16* - 5 • 4 ( - 9 • 8* + 3* • 4*) . (39) 
3d 



Considering the initial condition g\ = 0, we can solve 
Eq. (|39|) inductively leading to 



9t = 7^2* (158 + 112 • 16* - 270 ■ 8* + 210t • 4* 
945 v 



We now evaluate h t using a recursive method: 

2 f -l 

ht = 2* E ft(m)^(2*-m) 

m—l 
2*-l 



(40) 



2* 53 , 2 (m) 



m — l 

r.t-1 



534*+ 53 ^(fc) 



fe=i fc=i 

2 4 *- 1 + 2/i t _i , 



(41) 



where we have used the the symmetry v t (m) = v t (2* — m) . 
Since hi = 8, Eq. (|41 [) is solved inductively: 



(42) 



h t = 2 f+2 (8'-l)/7. 



To find an expression for A\' A , now the only thing left is to 
evaluate yt, which can be calculated using the same way 
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